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Compared  with  Mie  scattering  theory,  Airy  rainbow  theory  clearly  miscalculates  some  monochromatic 
details  of  scattering  by  small  water  drops.  Yet  when  monodisperse  Airy  theory  is  measured  by  percep¬ 
tual  (rather  than  purely  physical)  standards  such  as  chromaticity  and  luminance  contrast,  it  differs  very 
little  from  Mie  theory.  Considering  only  the  angular  positions  of  luminance  extrema,  Airy  theory’s 
errors  are  largest  for  small  droplets  such  as  those  that  dominate  cloudbows  and  fogbows.  However, 
integrating  over  a  realistic  drop-size  distribution  for  these  bows  eliminates  most  perceptible  color  and 
luminance  differences  between  the  two  theories.  ©  1998  Optical  Society  of  America 
OCIS  codes:  010.1290,  290.4020,  090.5640,  260.5430. 


1.  Introduction 

A  commonplace  of  optics  history  is  that  Newton’s 
(and  Descartes’)  geometrical  optics  dominated  18th- 
century  rainbow  theory.  Similarly,  George  Airy’s 
1838  interference/diffraction  theory1  built  on 
Thomas  Young’s  work  and  held  sway  for  several  de¬ 
cades  in  the  19th  century,  although  not  without  early 
detractors.  The  ascendance  of  these  theories  de¬ 
pended  on  their  ability  to  explain  naked-eye  features 
of  the  natural  rainbow  (i.e.,  bows  seen  in  rain,  clouds, 
fog,  or  spray).  However,  Airy  theory  was  soon  held 
to  a  different  standard:  predicting  the  angular  po¬ 
sitions  of  intensity  maxima  and  minima  for  spheres 
and  cylinders  illuminated  by  nearly  monochromatic 
light.2  Despite  early  success  with  such  predictions, 
Airy  theory  was  found  wanting  by  1888,  when  one 
experimenter  said  that  his  measurements  showed  it 
to  be  “but  a  first  approximation.”3 

Today  Airy  theory  is  seldom  compared  with  exact¬ 
ing  measurements,4  but  rather  with  the  intensity  dis¬ 
tribution  functions  of  Mie  theory.5  Here  Airy  theory 
may  not  position  the  rainbow  intensity  extrema  cor¬ 
rectly,  sometimes  turning  monochromatic  maxima 
into  minima  and  vice  versa.6’7  These  discrepancies 
are  most  noticeable  for  small  drop  sizes  and  parallel- 
polarized  rainbow  light.8  The  limitations  of  Airy’s 
cubic  wave-front  approximation  prompted  van  de 
Hulst  to  conclude  in  1957  that  the  “validity  of  Airy’s 
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theory  is  thus  limited  to  [size  parameters]  x  >  5000, 
or  with  light  of  X/2tt  =  0.1  p,  to  drops  with  radii  >¥z 
mm.”9  Naturally  van  de  Hulst  did  not  have  the  lux¬ 
ury  of  today’s  computer  resources,  so  checking  his 
analytic  claim  numerically  would  have  been  difficult. 
Yet,  even  in  1979,  Mobbs  cited  van  de  Hulst’s  claim  as 
one  reason  for  developing  his  own  rainbow  theory.10 

That  same  year,  however,  Konnen  and  de  Boer 
extended  Airy  theory  to  include  rainbow  polarization 
and  found  that  their  visible-wavelength  results  were 
reliable  at  drop  radii  as  small  as  —0.14  mm.11  Fur¬ 
thermore,  in  1977  Nussenzveig  described  monochro¬ 
matic  modeling  in  which  Airy’s  dominant  1 -polarized 
component  “requires  only  small  corrections  within 
the  primary  bow,  and  its  errors  become  appreciable 
only  in  the  region  of  the  supernumerary  arcs.”  Still, 
he  correctly  noted  that  “Airy’s  approximation  fails 
badly”  for  the  ||-polarized  primary.12  In  fairness, 
Wang  and  van  de  Hulst  have  recently  qualified  the 
earlier  proscription:  “We  found  that,  contrary  to 
what  has  often  been  thought ,  .  .  .  ,  Airy  theory  starts 
to  be  useful  at  relatively  small  sizes.  For  a  drop  size 
of  0.1  mm  it  already  represents  the  main  maxima  of 
the  primary  .  .  .  and  secondary  .  .  .  rainbows  quite 
well”  (emphasis  added).13 

The  italics  suggest  an  entrenched  conventional 
wisdom.  Certainly  the  implications  of  van  de 
Hulst’s  and  Nussenzveig’s  influential  earlier  analy¬ 
ses  are  clear:  (1)  Airy  theory  is  of  limited  use  in 
analyzing  the  rainbow,  and  (2)  Mie  theory  (or  at  least 
a  theory  more  sophisticated  than  Airy’s)  is  necessary 
for  quantitatively  reliable  rainbow  studies.  Kon¬ 
nen,  de  Boer,  Sassen,  Wang,  and  van  de  Hulst  have 
partially  exonerated  Airy  theory,  but  a  whiff  of  the 
unsavory  remains.  Thus  some  fresh  questions 
about  rainbow  verisimilitude  are  worth  asking. 
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Fig.  1.  Comparison  of  Mie  and  Airy  intensity  distribution  func¬ 
tions  at  size  parameter  x  =  628.3  for  wavelength  X  =  0.5  |xm  and 
a  drop  radius  of  r  =  50  |jum.  In  Figs.  1  and  2,  Crms  is  the  root- 
mean-square  distance  between  the  two  theories’  intensities  aver¬ 
aged  over  deviation  angle  0  =  137°-145°. 


First,  do  the  monochromatic  shortcomings  of  Airy 
theory  make  it  suspect  in  modeling  what  the  naked- 
eye  observer  sees  in  the  natural  rainbow?  In  other 
words,  should  Airy  theory  be  our  rainbow  theory  of 
last  resort?  Second,  how  perceptible  are  Airy  theo¬ 
ry's  photometric  and  colorimetric  departures  from 
Mie  theory? 

2.  Setting  the  Theoretical  Stage 

Even  a  cursory  comparison  of  Mie  and  Airy  theories 
reveals  distinct  differences  between  their  intensity 
distribution  functions  at  a  given  wavelength  and  drop 
size.  Figure  1  shows  this  difference  in  primary  rain¬ 
bows  at  wavelength  X  =  0.5  fxm  for  water  drops  with 
radii  r  =  50  p, m  (x  =  2'rrr/X  =  628.3).  Before  begin¬ 
ning  our  comparison,  however,  we  first  must  consider 
some  issues  of  terminology  and  scaling.  Figure  l's 
abscissa  poses  the  first  problem — Airy  theory  is  usu¬ 
ally  couched  in  terms  of  deviation  angle  0,  whereas 
Mie  theory  uses  scattering  angle  ©.  In  primary 
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deviation  angle  6 

(0  =  scattering  angle  0) 

Fig.  2.  Comparison  of  Mie  and  Airy  intensity  distribution  func¬ 
tions  at  x  ~  1885  for  X  =  0.5  |jim  and  r  =  150  fim. 


wavelength  (nm) 


Fig.  3.  Relative  solar  spectral  radiance  at  a  Sun  elevation  of 
—45°,  derived  from  measurements  at  University  Park,  Pa.,  on  5 
October  1987.  This  is  the  spectral  illuminant  assumed  in  calcu¬ 
lating  Figs.  4-25  and  27-30,  and  its  CIE  1976  uniform  chroma- 
ticity  scale  coordinates  are  u'  =  0.1986  and  v'  =  0.4713. 


bows,  0  and  ©  are  the  same,  but  in  the  secondaries  © 
=  360°  -  0.  For  consistency,  I  use  deviation  angle 
throughout  this  paper,  mindful  that  Mie  theory  does 
not  invoke  deviated  light  rays. 

The  ordinate  in  Fig.  1  also  requires  the  merging  of 
two  different  systems,  this  time  radiometric.  To  cal¬ 
culate  Airy  theory's  intensities,  I  use  the  formulas  of 
Humphreys,14  Tricker,15  and  Konnen  and  de  Boer.16 
Mie  theory  intensities  are  based  on  an  algorithm  by 
Bohren  and  Huffman  and  include  the  effects  of  ex¬ 
ternal  reflections.17  For  a  given  X,  r,  and  range  of  0, 
the  two  theories  can  yield  quite  different  intensity 
maxima,  so  I  must  normalize  one  theory's  results  in 
order  to  compare  them  with  the  other’s.  Here  I  scale 


deviation  angle  0 

(0  =  scattering  angle  0) 

Fig.  4.  Normalized  Mie  and  Airy  theory  primary  luminances  as 
functions  of  deviation  angle  0  for  a  single  500-fxm-radius  raindrop 
(both  polarizations).  The  two  theories’  intensity  distribution 
functions  are  convolved  with  Fig.  3’s  illuminant  to  produce  spec¬ 
trally  integrated  luminances.  These  luminances  are  plotted  in 
arbitrary  units  that  are  also  used  in  Figs.  5-7  and  18-21. 
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deviation  angle  0 

(0  =  scattering  angle  0) 

Fig.  5.  Normalized  Mie  and  Airy  theory  primary  luminances  as 
functions  of  deviation  angle  0  for  a  single  150-fxm-radius  drizzle 
drop  (both  polarizations). 


Airy  theory  curves  by  the  normalizing  factor  /max, Mie/ 
4iax,Airy>  where  7max> Mie  and  7max> Airy  are  the  respec- 
tive  intensity  (or  luminance  Lv )  maxima  for  the  two 
Mie  and  Airy  theory  curves  being  compared.  Thus 
in  every  figure  that  compares  intensity  (or  lumi¬ 
nance),  my  Airy  and  Mie  maxima  are  the  same. 
Other  normalizations  are  possible,  but  they  do  not 
alter  my  qualitative  conclusions. 

In  Fig.  1,  the  Mie  ripple  structure  (both  polariza¬ 
tions)  is  evident  as  it  oscillates  around  the  smoother 
Airy  intensity  curve  (_L  polarization  only).  Such 
monochromatic  comparisons  usually  illustrate  the 
failings  of  Airy  theory,  especially  of  its  ||-polarized 
component.18  However,  Fig.  1  also  shows  how  the 
Airy  J_ -polarized  component  approximately  follows 
the  Mie  extrema.  At  each  0  and  corresponding  ©  in 
Fig.  1, 1  quantify  this  difference  as  an  Airy-Mie  con- 
trast  C,  where  C  =  (7^  -  7Mie)/7Mie.  In  Fig.  1,  the 
root-mean-square  (rms)  contrast  difference  between 
the  two  theories  is  0.4186.  Contrary  to  expecta- 


deviation  angle  0 

(0  =  scattering  angle  0) 


Fig.  6.  Normalized  Mie  and  Airy  theory  primary  luminances  as 
functions  of  deviation  angle  0  for  a  single  50-|jim-radius  cloud  drop 
(both  polarizations). 


deviation  angle  0 

(0  =  scattering  angle  0) 

Fig.  7.  Normalized  Mie  and  Airy  theory  primary  luminances  as  a 
function  of  deviation  angle  0  for  a  single  10-|jim-radius  cloud  drop 
(both  polarizations). 


tions,  this  rms  difference  rises  to  0.5061  for  Fig.  2’s 
larger  drop  radius  (150  |xm),  an  increase  attributable 
both  to  larger  excursions  in  the  Mie  ripple  structure19 
and  smaller  Airy  minima.  Cast  in  these  monochro¬ 
matic  terms,  Airy  theory  is  indeed  a  poor  second  to 
Mie  theory. 

However,  if  we  compare  the  two  theories’  colors  and 
luminances,  the  differences  are  subtler.  My  com¬ 
parisons  make  some  assumptions  that  are  conserva¬ 
tive  (i.e.,  they  preserve  some  Mie  scattering  details) 
and  sometimes  literally  unnatural  (i.e.,  they  do  not 
include  all  factors  affecting  natural  bows).  Thus  I 
am  not  relentlessly  smoothing  away  Mie  theory  de¬ 
tails,  but  instead  examining  whether  they  result  in 


Fig.  8.  Portion  of  the  CIE  1976  UCS  diagram,  showing  the  chro- 
maticity  of  Fig.  3’s  illuminant  (+)  and  separate  77'(0)>  i/(0)  chro- 
maticity  curves  for  Airy  and  Mie  theory  primary  cloudbows  at 
10-|xm  radius  (both  polarizations).  The  two  curves’  mean  colori¬ 
metric  separation  A {u' ,  v')  =  0.0080,  compared  with  a  mean  Mac- 
Adam  JND  of  0.004478  in  the  u' ,  v'  region  spanned  by  the  Mie  and 
the  Airy  primaries. 
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Fig.  9.  Close-up  view  of  Fig.  8.  Note  that  the  u',  v'  scaling  is 
isotropic  in  this  and  all  subsequent  UCS  diagrams. 


Fig.  11.  Chromaticity  curves  for  the  50-|jini-radius  Airy  and  Mie 
theory  primaries  differ  by  A (u\  v')  =  0.00571  (both  polarizations). 


rainbow  features  visible  to  naked-eye  observers.  In 
all  that  follows,  bear  in  mind  that  I  am  not  empha¬ 
sizing  the  two  theories'  electromagnetic  or  mathe¬ 
matical  details,  but  instead  their  visible  differences. 

First,  unlike  Figs.  1  and  2,  all  subsequent  figures 
describe  the  convolution  of  a  rainbow  intensity  dis¬ 
tribution  function  and  a  particular  illuminant.  Fig¬ 
ure  3  shows  this  illuminant’s  spectrum,  which 
corresponds  to  sunlight  measured  at  the  Earth’s  sur¬ 
face  when  the  Sun’s  elevation  is  —45°.  Obviously 


u" 

Fig.  10.  Perspective  view  that  combines  chromaticities  (Fig.  9) 
and  luminances  (Fig.  7)  for  Mie  and  Airy  theory  10-|jim-radius 
primary  cloudbows  (both  polarizations).  The  combined  lumi¬ 
nance  and  chromaticity  curves  are  also  shown  projected  on  the  u' , 
v'  chromaticity  plane. 


very  few  rainbows  will  be  seen  for  a  45°  Sun  eleva¬ 
tion,20  but  my  point  here  is  simply  to  choose  a  natural 
illuminant  that  is  not  highly  chromatic.  Second,  all 
subsequent  figures  include  the  smoothing  effects  of 
sunlight’s  approximately  0.5°  angular  divergence.21 
Third,  all  comparisons  include  both  the  4-  and  the 
||-polarized  components  of  rainbow  light,  just  as  the 
naked-eye  observer  must.  However,  except  in  two 
cases,  I  do  make  the  unnatural  assumption  that  only 
a  single,  spherical  drop  generates  the  bows.  In  other 
words,  because  I  integrate  over  a  drop-size  distribu¬ 
tion  only  rarely,  my  canonical  comparison  is  between 
monodisperse  luminances  and  chromaticities  for  Mie 
and  Airy  theories. 


Fig.  12.  Chromaticity  curves  for  the  150-fxm-radius  Airy  and  Mie 
theory  primaries  differ  by  A (u' ,  v')  =  0.01037  (both  polarizations). 
Breaks  in  the  Airy  curve  mark  the  positions  of  its  primary  and  first 
supernumerary. 
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Fig.  13.  Chromaticity  curves  for  the  500-fxm-radius  Airy  and  Mie 
theory  primaries  differ  by  A (u' ,  v')  =  0.01280  (both  polarizations). 

3.  Monodisperse  Comparisons  of  Mie  and  Airy 
Primaries 


and  10  fxm.22  For  a  large  500-fxm-radius  drop,  Fig. 
4  shows  the  expected  close  agreement  between  Airy 
and  Mie  theories.23  Now  define  luminance  contrast 
as  Cv  =  [(Lt,  Ajry  -  LvMie)/Lv  Mie],  At  a  given  drop 
radius,  Cv  is  averaged  over  N  discrete  Ly(0j)  to  yield 
the  rms  contrast  Crms,  where 


n 

^rms 


N 

2 

i=l 


^v,Airy(^i)  ^v,Mi  e(®i) 

^v,Mi  eWJ 


2  'j  1/2 

/N 


(1) 


and  N  =  383.  In  Fig.  4,  Crms  =  0.1327,  a  much 
smaller  value  than  those  for  the  monochromatic  Figs. 
1  and  2.  Not  surprisingly,  integration  over  the  vis¬ 
ible  spectrum  eliminates  much  of  the  Airy-Mie  lumi¬ 
nance  difference.  Like  most  contrast  measures,  Cv 
scales  luminance  differences  in  terms  of  some  refer¬ 
ence  luminance,  which  here  is  Lv  Mie.  Thus  Crms 
shows  relative  departures  from  our  Mie  rainbow 
standard,  and  those  departures  are  averaged  across 
the  deviation  angle. 

Note  that  integration  over  the  Sun's  diameter  in 
Fig.  4  reduces  the  original  L„(0)  to  the  range  137.23° 
<  0  <  144.76°.  This  integration  takes  the  form 


We  begin  by  comparing  the  spectrally  integrated  rel¬ 
ative  luminances  L„(0)  of  single-droplet  primary  rain¬ 
bows  and  cloudbows  with  drop  radii  of  500,  150,  50, 


LM  = 


fQ+w 

Q-ftV 

1_  - 

Jd-w 

\  w  )  J 

1/2 


L,(fi)dfi,  (2) 


Fig.  14.  Map  of  Airy  theory  colors  for  monodisperse  primary  rainbows  and  cloudbows.  The  map  assumes  (1)  Fig.  3’s  solar  spectrum  as 
the  illuminant,  (2)  Eq.  (2)’s  sun-width  smoothing  filter,  (3)  both  rainbow  polarizations,  (4)  spherical,  nonabsorbing  water  drops.  At  each 
drop  radius,  colors’  luminances  are  normalized  by  the  maximum  luminance  for  that  radius.  A  colorimetrically  calibrated  version  of  this 
figure  can  be  seen  at  the  www  address  given  in  the  acknowledgment. 
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Fig.  15.  Map  of  Mie  theory  colors  for  monodisperse  primary  rainbows  and  cloudbows.  Figure  14’s  assumptions  (l)-(4)  also  hold  here. 
A  colorimetrically  calibrated  version  of  this  figure  can  be  seen  at  the  www  address  given  in  the  acknowledgment. 


where  LS(Q)  is  the  smoothed  luminance  and  w  is  the 
solar  radius  of  —0.25°  [all  Eq.  (2)  angles  are  in  radi¬ 
ans].  Equation  (2)’s  square-root  term  is  propor¬ 
tional  to  the  Sun’s  angular  width  at  each  radial  angle 
fl  across  the  Sun’s  disk.  In  effect,  at  each  0  we  ap¬ 
proximate  the  finite-width  sun  as  the  sum  of  many 
point-source  suns,  each  of  which  is  weighted  by  the 
real  Sun’s  angular  width  at  radial  angle  fi.  In  turn, 
each  point-source  sun  contributes  its  own  rainbow 
Lv(ft)  to  the  observed  Ls( 0).  Integrals  similar  to  Eq. 
(2)  smooth  the  Mie  and  the  Airy  chromaticities  dis¬ 
cussed  below. 

How  significant  visually  are  Fig.  4’s  Crms?  Small 
tick  marks  shown  in  the  ±  threshold  contrast  boxes  of 
Figs.  4-7  provide  one  indication.  Each  tick  mark 
spans  a  luminance  range  that  is  ±2%  of  its  mean 
value.  If  a  ±2%  luminance  change  is  a  just- 
noticeable  difference24  (JND),  then  Fig.  4’s  tick  mark 
is  a  graphical  measure  of  threshold  contrast.  Note 
that  Fig.  4’s  tick  mark  will  have  the  same  length 
anywhere  along  its  logarithmic  ordinate,  although 
that  length  will  change  when  the  ordinate’s  range 
changes.  Gauged  in  terms  of  contrast,  then,  most 
Airy-Mie  luminance  differences  in  Fig.  4  are  sub¬ 
threshold  (i.e.,  they  are  invisible).  Even  in  the  su¬ 
pernumeraries  (0  =  141°-143.5°),  Fig.  4’s  Airy  theory 
errors  usually  are  invisible.  The  largest  disagree¬ 
ment  occurs  in  Alexander’s  dark  band  (0  <  137.6°), 


where  Mie  theory  consistently  predicts  more  light.25 
Yet  the  overall  close  agreement  in  Fig.  4  is  not  sur¬ 
prising,  because  the  mean  visible-wavelength  size  pa¬ 
rameter  x  =  6004  here,  which  exceeds  van  de  Hulst’s 
1957  lower  limit  on  Airy  theory’s  validity. 

Figure  5  shows  Airy  and  Mie  luminances  for  a 
150-|xm  drizzle  drop.  Crms  has  decreased  slightly  to 
0.1122,  as  has  the  luminance  difference  in  Alex¬ 
ander’s  dark  band.  However,  the  Airy  underesti¬ 
mates  of  luminance  in  the  supernumerary  minima 
now  are  near  or  above  threshold.  Supernumerary 
maxima  are  still  indistinguishably  different  in  the 
two  theories.  Although  x  (150  fxm)  =  1801,  Airy 
departures  from  Mie  theory  still  are  not  pervasive. 
However,  in  Fig.  6’s  50-fxm  cloudbow  (x  =  600),  small 
angular  offsets  of  the  Airy  supernumeraries  from 
their  Mie  counterparts  are  now  evident  even  in  the 
first  supernumerary.  Airy  supernumerary  minima 
are  now  noticeably  darker  than  those  for  Mie  theory, 
as  the  increase  of  Crms  to  0.1240  suggests. 

For  a  small  cloud  drop  (10-|xm  radius,  x  =  120),  the 
two  bows  have  broadened  sufficiently  in  Fig.  7  that  no 
supernumeraries  are  visible,  although  residual 
ripple-structure  fluctuations  are  still  evident  in  the 
Mie  cloudbow.  (These  luminance  ripples  do  have  a 
color  counterpart,  as  shown  below.)  Although  Crms 
is  smaller  (0.1082)  than  that  for  Fig.  6’s  larger  cloud 
drop,  Fig.  7’s  Mie  and  Airy  curves  are  even  less  con- 
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Fig.  16.  Map  of  Airy  theory  colors  for  monodisperse  secondary  rainbows  and  cloudbows.  Figure  14’s  assumptions  (l)-(4)  also  hold  here. 
A  colorimetrically  calibrated  version  of  this  figure  can  be  seen  at  the  www  address  given  in  the  acknowledgment. 


gruent.  This  is  largely  due  to  the  Airy  cloudbow’s 
shift  to  larger  0,  which  places  its  peak  luminance  —1° 
closer  to  the  antisolar  point  than  for  the  Mie  10-|xm 
cloudbow.  Although  this  shift  is  obvious  graphi¬ 
cally,  observing  it  in  real  cloudbows  is  problematic,  as 
shown  below. 

Now  we  consider  chromaticity  differences  between 
Mie  and  Airy  theories  for  the  same  four  monodisperse 
primaries  (r  =  10,  50,  150,  and  500  p,m).  All  chro- 
maticities  are  calculated  by  summations  from  380  to 
700  nm  in  5-nm  steps.26  Figure  8  shows  a  portion  of 
the  CIE  1976  uniform  chromaticity  scale  (UCS)  dia¬ 
gram.  In  it,  the  illuminant’s  color  is  marked  with  a 
+ ;  two  chromaticity  curves  trace  the  Mie  and  the  Airy 
u' ,  v'  across  0  for  r  =  10  p,m.  One  measure  of  color 
difference  between  the  two  theories  is  the  colorimet¬ 
ric  distance  A (ur,  v').  At  any  deviation  angle  0, 

A (w\  v')  =  {[M'(6)Mie  -  «'(6)Airy]2  +  [l>'(0)Mie 

-  u'(e)Airy]2}1/2,  (3) 

and  the  average  A (ur,  v')  is  calculated  over  the  0 
range  being  considered.  In  Fig.  8,  A(z/',  vf)  =  0.008, 
—75%  greater  than  the  mean  MacAdam  JND27  of 
0.004478  in  the  UCS  region  spanned  by  the  Mie  and 
the  Airy  primaries.  Like  the  contrast  tick  mark 
used  in  Figs.  4-7,  the  MacAdam  JND  drawn  in  Fig. 
8  serves  as  a  graphical  ruler  of  threshold  difference. 


Thus  for  a  10-|xm-radius  drop,  Mie  and  Airy  theo¬ 
ries  produce  nearly  achromatic  primaries  that  are 
colorimetrically  distinguishable,  at  least  on  their 
exteriors.  Figure  9’s  close-up  view  of  Fig.  8  shows 
the  differences  between  Mie  theory’s  chromaticity 
loops  and  the  simpler  hook-shaped  curve  of  Airy 
theory.  Figure  10  joins  Fig.  7’s  luminances  and 
Fig.  9’s  chromaticities  to  give  a  perspective  view  of 
their  combined  variation.  Certainly  the  Mie  loops 
are  physically  justifiable,  yet  their  complexity 
seems  at  odds  with  the  simple  color  transitions  seen 
in  real  cloudbows.28 

Similar  loops  or  wiggles  appear  on  the  exterior  of 
Fig.  ll’s  Mie  chromaticity  curve  for  a  50- [xm  cloud 
drop.  However,  here  the  wiggles  are  much  less 
prominent.  In  fact,  now  the  Mie  and  Airy  chroma¬ 
ticities  correspond  much  better,  as  is  evident  both 
graphically  and  numerically  [A (ur,  v')  =  0.00571]. 
Only  occasionally  do  the  two  theories  disagree  by 
more  than  1  JND.  Figure  12  shows  that  at  150 
[xm,  A (V,  v')  has  increased  to  0.01037,  contrary  to 
the  conventional  wisdom  that  differences  between 
Airy  and  Mie  theory  must  always  decrease  with 
increasing  drop  size.29  Another  plausible  assump¬ 
tion  upset  by  Fig.  12  is  that  rainbow  luminance  and 
chromaticity  extremes  necessarily  coincide.  The 
Airy  primary  and  first  supernumerary  (i.e.,  the  two 
leftmost  maxima  in  Fig.  5)  are  indicated  by  small 
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Fig.  17.  Map  of  Mie  theory  colors  for  monodisperse  secondary  rainbows  and  cloudbows.  Figure  14’s  assumptions  (l)-(4)  also  hold  here. 
A  colorimetrically  calibrated  version  of  this  figure  can  be  seen  at  the  www  address  given  in  the  acknowledgment. 


breaks  in  Fig.  12’s  chromaticity  curve.  Clearly 
these  two  luminance  maxima  do  not  have  purer 
colors  than  their  neighbors.  In  fact,  Airy  theory 
predicts  that  the  primary’s  purest  colors  will  be 
very  dark  reds  on  its  exterior.  Mie  theory  makes 
these  colors  both  slightly  brighter  (see  Fig.  5)  and 
less  pure  than  those  slightly  closer  to  the  primary 
maximum. 

The  divergence  between  Mie  and  Airy  colors  out¬ 
side  the  primary  is  even  greater  at  a  5Q0-|xm  radius 
(Fig.  13),  and  it  accounts  for  most  of  the  A (u',  v’) 
reported  there  (0.01280).  In  fact,  eliminating  the 
smallest  deviation  angles  (0  <  137.95°)  from  A (u' ,  v') 
reduces  it  threefold  to  0.00373.  Although  the  Mie 
theory  desaturation  outside  Fig.  13’s  primary  is  quite 
dramatic,  we  are  unlikely  to  see  it  in  nature,  as  rain¬ 
bow  light  will  be  additively  mixed  with  background 
light.  Because  luminance  on  the  Mie  500- [xm  pri¬ 
mary’s  exterior  is  —63  times  less  than  at  its  peak 
(Fig.  4),  background  light  from  clouds  will  almost 
always  dominate  colors  outside  the  natural  bow.  Fi¬ 
nally,  note  that  the  purity  of  Fig.  13’s  Airy  primary 
has  increased  dramatically  from  that  of  Fig.  12.  Al¬ 
though  the  dominant  wavelength  is  —566  nm  for 
both,  the  Airy  500-fxm  primary  is  both  yellower 
(82.3%  versus  34.7%  purity)  and  brighter  than  its 
150-|xm  counterpart. 


4.  Color  Maps  of  Mie  and  Airy  Primaries 

As  useful  as  Figs.  4-13  are  in  plotting  selected  dif¬ 
ferences  between  the  Mie  and  the  Airy  primaries, 
they  do  not  convey  any  overall  visual  sense  of  these 
predicted  rainbows.  Because  our  benchmark  here  is 
the  two  theories’  perceptual  disagreement,  we  also 
need  a  purely  visual  comparison.  Figures  14  and  15 
provide  it  in  the  form  of  color  maps  of  the  Airy  and 
Mie  primaries,  respectively;  Figs.  16  and  17  map  the 
Airy  and  Mie  secondaries.  The  droplet  radius  in¬ 
creases  logarithmically  along  each  map’s  abscissa, 
whereas  the  deviation  angle  varies  linearly  along  the 
ordinate.  All  maps  are  arranged  so  that  their  top- 
to-bottom  color  sequence  is  the  same  as  that  seen  at 
the  summit  of  the  natural  primary  or  secondary. 

To  make  each  map,  standard  projective  geometry 
techniques  are  used  to  convert  rainbow  chromatici- 
ties  to  their  red- green- blue  equivalents  on  a  com¬ 
puter’s  calibrated  color  monitor.30  Figures  14-17 
are  mapped  so  that  their  white  corresponds  to  Fig.  3’s 
achromatic  u' ,  v' .  (Your  achromatic  u\  v'  depends 
on  the  illuminant  with  which  you  view  Figs.  14-17.) 
Neither  a  computer  monitor  nor  the  printed  page  can 
adequately  reproduce  colors  throughout  the  lumi¬ 
nance  dynamic  range  shown  in  Figs.  4-7.  Thus 
each  column  of  monodisperse  colors  in  Figs.  14-17  is 
normalized  by  the  maximum  luminance  found  at  that 
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drop  size.  In  other  words,  Figs.  14-17  do  not  show 
the  enormous  range  of  rainbow  luminances  evident 
in  Figs.  4-7.  Like  Figs.  4-13,  Figs.  14-17  assume 

(1)  Fig.  3’s  solar  spectrum  as  the  illuminant,  (2)  Eq. 

(2) ’s  sun-width  smoothing  filter,  (3)  both  rainbow  po¬ 
larizations,  (4)  spherical,  nonabsorbing  water  drops, 
and  (5)  monodisperse  colors.  A  monodisperse  com¬ 
parison  is,  of  course,  far  stricter  than  that  possible  in 
the  natural  bow. 

The  most  striking  feature  of  Figs.  14  and  15  is  their 
essential  similarity.  Even  at  cloud-drop  sizes,  the 
positions  of  the  Airy  and  the  Mie  primary  maxima 
appear  nearly  identical  (but  see  Fig.  7).  In  fact,  the 
only  literally  unnatural  feature  seems  to  be  the  Mie 
map’s  ripple  structure,  which  appears  as  a  subtle 
color  marbling  on  the  cloudbow  primaries  and  their 
supernumeraries.  As  drop  radius  increases,  in  both 
Figs.  14  and  15  the  primary  and  the  supernumeraries 
converge  into  ever-narrower  bands.  At  drop  radii 
>0.4  mm,  Mie  and  Airy  supernumeraries  are  nearly 
achromatic  and  individual  supernumeraries  are  of¬ 
ten  blurred  because  of  Eq.  (2)’s  smoothing.  At  those 
same  sizes,  we  can  also  see  that  Alexander’s  dark 
band  is  brighter  for  Mie  than  Airy  theory  (see  Fig.  4). 

Quantitatively,  a  pixel-by-pixel  comparison  of  Figs. 
14  and  15  shows  that  only  5.32%  of  all  Mie  and  Airy 
rainbow  colors  have  A (u',  v')  >  0.02.31  Color  dis¬ 
tances  >0.02  are  largely  confined  to  Alexander’s  dark 
band  at  large  drop  sizes,  where  very  low  luminances 
render  these  large  mismatches  essentially  meaning- 
less.  For  pixels  with  A(^',  vr)  <  0.02,  A (u* ,  v')  = 
0.005885,  and  its  standard  deviation  =  0.003873. 
Thus  the  average  Mie-Airy  color  distance  only 
slightly  exceeds  the  primary  rainbow  JND.  In  fact, 
over  41%  of  all  Airy  colors  here  are  within  one  JND  of 
their  Mie  counterparts  (i.e.,  compare  Airy  and  Mie 
colors  at  the  same  0  and  r).  For  the  difficult  condi¬ 
tions  of  outdoor  color  matching,  even  2  JND’s  are  a 
defensible  threshold,  in  which  case  over  79%  of  Airy 
theory  colors  would  be  confused  with  their  Mie  coun¬ 
terparts.  Any  presumed  colorimetric  chasm  be¬ 
tween  the  two  theories  now  looks  considerably 
smaller.  In  fact,  much  of  the  Airy-Mie  A (u’,  v')  is 
due  to  the  Mie  ripple  structure,  which  is  not  seen  in 
natural  bows. 

A  similar  comparison  of  luminance  contrast  in 
Figs.  14  and  15  yields  a  mean  Cv  of  -0.02267  with  a 
standard  deviation  of  0.1141.  For  a  threshold  con¬ 
trast  \CV\  =  0.02,  some  31%  of  the  two  maps’  pixels 
would  be  confused.  Furthermore,  the  positions  of 
primary  and  supernumerary  luminance  peaks  and 
valleys  in  Figs.  14  and  15  are  quite  close  at  most  drop 
sizes.  This  means  that  subthreshold  contrasts  occur 
at  all  drop  sizes,  not  just  at  raindrop  sizes.  Quanti¬ 
tatively,  monodisperse  Airy  primaries  often  do  depart 
perceptibly  from  their  Mie  counterparts.  Yet  quali¬ 
tatively,  these  Airy  primaries  look  more  natural  than 
do  the  Mie  rainbows  and  cloudbows.  By  this  I  mean 
that  each  monodisperse  Airy  rainbow  more  closely 
resembles  a  naturally  occurring  polydisperse  rain¬ 
bow  or  cloudbow  than  does  its  Mie  counterpart. 


229.5°  231.5°  233.5°  235.5°  237.5° 

deviation  angle  0 

(scattering  angle  ©  =  360°  -  0) 

Fig.  18.  Normalized  Mie  and  Airy  theory  secondary  luminances 
as  functions  of  deviation  angle  0  for  a  single  500-|jim-radius  rain¬ 
drop  (both  polarizations).  The  illuminant  is  Fig.  3’s  sunlight  spec¬ 
trum. 


5.  Monodisperse  Comparisons  of  Mie  and  Airy 
Secondaries 

Mie  and  Airy  secondaries  might  be  expected  to  differ 
more  than  the  primaries,  and  Figs.  18-21  do  indeed 
show  greater  luminance  differences  in  the  secondar¬ 
ies.  Aside  from  a  different  range  of  0,  all  other  Mie 
and  Airy  theory  conditions  are  unchanged  from  those 
of  Figs.  4-7.32  In  the  secondary,  note  that  scatter¬ 
ing  angle  ©  =  360°  -  0  and  that  0  itself  is  now  located 
0  -  180°  from  the  antisolar  point.  In  Fig.  18,  Crms 
has  more  than  doubled  to  0.3108  from  Fig.  4’s  value, 
with  most  of  the  increase  due  to  Airy  underestimates 
of  luminances  in  Alexander’s  dark  band  (0  <  231°). 
Similar  contrast  increases  are  evident  in  Figs.  19-21, 
where  Crms  ranges  from  0.2839  to  0.2946.  As  these 
figures’  tick  marks  indicate,  contrast  between  the 
monodisperse  Mie  and  Airy  secondaries  usually  ex¬ 
ceeds  threshold.  Next  we  examine  the  angular 
shifts  of  Airy  maxima  and  minima  in  Figs.  19-21 


deviation  angle  0 

(scattering  angle  0  =  360°  -  0) 


Fig.  19.  Normalized  Mie  and  Airy  theory  secondary  luminances 
as  functions  of  deviation  angle  0  for  a  single  150-|jim-radius  drizzle 
drop  (both  polarizations). 
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Fig.  20.  Normalized  Mie  and  Airy  theory  secondary  luminances 
as  functions  of  deviation  angle  0  for  a  single  50-|jim-radius  cloud 
drop  (both  polarizations). 


from  the  Mie  positions.  These  shifts  are  larger  than 
similar  Airy  displacements  in  the  primaries  (Figs. 
5-7). 

Comparing  Figs.  4  and  18,  note  that  the  MieiAiry 
luminance  ratio  in  Alexander's  dark  band  is  much 
larger  in  the  secondary  than  in  the  primary.  The 
same  is  true  at  other  drop  sizes,  meaning  that  Airy 
theory  consistently  makes  this  part  of  the  second¬ 
ary  too  dark.  Similarly,  the  MieiAiry  luminance 
ratio  is  slightly  larger  for  supernumerary  minima 
in  Fig.  19  than  in  Fig.  5.  (Also  note  that  the  Airy 
supernumerary  maxima  are  too  bright  in  Fig.  19.) 
Mathematically,  this  occurs  because  the  Airy  sec¬ 
ondaries'  _L-  and  ||-polarized  intensities  remain  in 
phase,  unlike  their  primary  counterparts.33  In 
other  words,  Airy  theory's  ||-polarized  component 
increases  minima  less  in  its  secondaries  than  in  its 
primaries,  where  the  ||-  and  the  T-polarized  compo¬ 
nents  are  almost  completely  out  of  phase.  Thus  a 
deficiency  of  the  Airy  ||-polarized  primary  (its  out- 


deviation  angle  9 

(scattering  angle  0  =  360°  -  0) 


Fig.  21.  Normalized  Mie  and  Airy  theory  secondary  luminances 
as  functions  of  deviation  angle  0  for  a  single  10-|jim-radius  cloud 
drop  (both  polarizations). 


u' 


Fig.  22.  Chromaticity  curves  for  the  10-fjim-radius  Airy  and  Mie 
theory  secondaries  differ  by  A (u' ,  v')  =  0.01585  (both  polariza¬ 
tions).  The  mean  MacAdam  JND  in  the  u' ,v'  region  spanned  by 
the  Mie  and  the  Airy  secondaries  =0.004773. 

of-phase  behavior)34  actually  improves  the  Airy  pri¬ 
mary’s  fit  to  Mie  theory  at  some  drop  sizes,  provided 
that  we  consider  both  polarizations  (e.g.,  Fig.  5). 

The  Mie  and  the  Airy  secondary  chromaticities  also 
yield  some  counterintuitive  results.  These  are 
shown  in  Figs.  22-25,  where  A (uf,  v’)  consistently 
increases  with  drop  size.  The  secondaries’  A (uf,  v') 
are  much  larger  than  the  primaries',  being  some  3  to 
8  times  greater  than  the  mean  MacAdam  JND  of 
0.004773  for  Mie  and  Airy  secondaries.  Yet  despite 
the  mathematical  and  physical  rigor  that  underlies 
the  Mie  chromaticity  curves,  they  appear  quite  liter¬ 
ally  unnatural.  In  Fig.  22,  the  10-fxm-radius  Mie 
chromaticities  trace  out  a  convoluted  curve  that  dif¬ 
fers  starkly  from  the  simple  arc  of  the  Airy  secondary 


Fig.  23.  Chromaticity  curves  for  the  50-|jim-radius  Airy  and  Mie 
theory  secondaries  differ  by  A (u' ,  v')  =  0.01944  (both  polariza¬ 
tions). 
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Fig.  24.  Chromaticity  curves  for  the  150-|xm-radius  Airy  and  Mie 
theory  secondaries  differ  by  A (u' ,  v')  =  0.03324  (both  polariza¬ 
tions). 


cloudbow.  Of  course,  at  this  drop  size  not  all  the 
secondary  is  visible  for  229.7°  <  0  <  237.3°,  but  the 
difference  between  the  two  theories  is  clear  enough. 

At  50-fxm-radius  (Fig.  23),  the  Airy  colors  resemble 
a  curved  frame  for  the  irregular  helix  of  Mie  chroma- 
ticities.  Near  the  center  of  the  Airy  curve,  color  dis¬ 
tances  are  subthreshold  more  often  than  at  10  fxm, 
yet  the  50-fxm  A (u' ,  v’)  is  greater  because  Airy  theory 
predicts  a  larger  color  gamut  (i.e.,  the  Airy  curve 
extends  beyond  the  Mie  curve).35  For  a  150-|xm- 
radius  drizzle  drop,  Fig.  24’s  Mie  and  Airy  chroma- 
ticities  at  first  look  similar  to  their  Fig.  12 
counterparts.  However,  A (u',  u')  is  3  times  larger  in 
Fig.  24,  a  consequence  of  the  Mie  secondary  tracing  a 
tortuous  path  within  its  smaller  gamut.  In  fact, 
these  Mie  wiggles  are  colorimetric  manifestations  of 


Fig.  17’s  skewed  color  bands  at  150  [xm,  themselves  a 
result  of  the  color  ripple  structure.  At  500-fxm  ra¬ 
dius  (Fig.  25),  the  Mie  secondary  chromaticities  rip¬ 
ple  only  occasionally,  and  like  Fig.  13’s  Mie 
primaries,  they  rapidly  desaturate  within  Alex¬ 
ander’s  dark  band.  Although  Fig.  25’s  Mie  and  Airy 
curves  are  hardly  congruent,  their  similar  shapes 
seem  as  noteworthy  as  their  large  A(z/',  v'). _ 

Given  the  secondaries’  dramatically  larger  A (u' ,  vf), 
Figs.  16  and  17  are  surprising.  Although  these  color 
maps  do  differ,  they  are  not  the  distant  relatives  that 
Figs.  18-25  suggest.  Beginning  at  r  ~  0.07  mm  in 
Fig.  16,  Airy  theory’s  underestimates  of  minimum  su¬ 
pernumerary  luminances  are  evident  (i.e.,  correspond¬ 
ing  luminance  minima  are  darker  than  in  Fig.  17). 
Also  evident  is  the  consistently  brighter  Alexander’s 
dark  band  in  Fig.  17’s  lower  right-hand  corner.  Colo- 
rimetrically,  12.1%  of  Fig.  16’s  Airy  pixels  are  within  1 
JND  of  corresponding  Mie  pixels  in  Fig.  17;  32.4%  are 
within  2  JND’s.  As  in  the  primary,  many  of  the  larg¬ 
est  Airy  colorimetric  errors  in  the  secondary  occur  at 
low  luminances  in  Alexander’s  dark  band,  where  they 
are  unlikely  to  be  visible.  Similarly,  significant  num¬ 
bers  of  subthreshold  A (uf,  v’)  and  Cv  are  found  at  all 
drop  sizes  in  Fig.  16,  although  they  do  occur  more 
frequently  at  larger  sizes. 

Although  we  know  objectively  that  some  Airy  su¬ 
pernumeraries  are  shifted  (e.g.,  Figs.  20-21),  gross 
discrepancies  are  difficult  to  see  in  Fig.  16.  In  part 
this  is  due  to  the  distracting  color  marbling  of  Fig. 
17’s  Mie  ripple  structure  below  r  ~  0.1  mm.  [Note 
that  the  disappearance  of  this  structure  at  larger 
radii  is  just  an  artifact  of  Fig.  17’s  (and  Fig.  15’s) 
angular  resolution.]  In  fact,  Fig.  17’s  color  marbling 
nearly  obscures  the  smaller  drops’  subtle  luminance 
peaks.  At  drizzle  and  raindrop  sizes  (r  >  0.1  mm), 
the  two  maps  show  that  Mie  and  Airy  secondaries  are 
almost  indistinguishable  visually.  Furthermore, 
despite  its  larger  quantitative  errors,  Fig.  16’s  Airy 


Fig.  25.  Chromaticity  curves  for  the  500-fjum-radius  Airy  and  Mie 
theory  secondaries  differ  by  A (u' ,  v')  =  0.04019  (both  polariza¬ 
tions). 


Fig.  26.  Deirmendjian  modified  gamma  distribution  for  drop 
sizes  typical  of  stratocumulus  or  fog,  l-40-|jim  radii.  This  poly¬ 
dispersion  produces  the  smoothed  Mie  and  Airy  theory  data  shown 
in  Figs.  27-30. 
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deviation  angle  0 

(0  =  scattering  angle  0) 


Fig.  27.  Normalized  luminances  for  Mie  and  Airy  primary  fog- 
bows  as  functions  of  deviation  angle  0  for  Fig.  26’s  drop-size  dis¬ 
tribution  (both  polarizations).  In  this  and  subsequent  figures,  the 
illuminant  is  Fig.  3’s  sunlight  spectrum. 


map  arguably  offers  a  clearer  picture  of  natural  sec¬ 
ondary  bows  than  does  Fig.  17.36  Like  its  primary 
counterpart  (Fig.  14) ,  Fig.  16  seems  the  more  useful 
qualitative  guide  to  rainbows  and  cloudbows  viewed 
with  the  naked  eye. 

6.  Toward  Real  Fogbows:  Mie  and  Airy 
Polydispersions 

As  instructive  as  Figs.  4-25  are,  they  still  lack  an 
important  natural  detail — smoothing  over  a  realistic 
drop-size  spectrum.  Because  Airy  theory  suppos¬ 
edly  differs  most  from  Mie  theory  at  small  drop  radii, 
I  integrate  over  a  poly  dispersion  of  cloud  droplets. 
Here  I  use  a  Deirmendjian  modified  gamma  distribu¬ 
tion  that  is  representative  of  drop-size  spectra  ob¬ 
served  in  stratocumulus  and  fog  (Fig.  26).37  Not 
surprisingly,  when  Mie  and  Airy  luminances  are 
weighted  by  Fig.  26’s  droplet  number  densities,  the 
resulting  primary  and  secondary  fogbows  resemble 
each  other  more  closely  than  do  the  monodisperse 


deviation  angle  0 

(scattering  angle  0  =  360°  -  0) 


Fig.  28.  Normalized  luminances  for  Mie  and  Airy  secondary  fog¬ 
bows  as  functions  of  deviation  angle  0  for  Fig.  26’s  drop-size  dis¬ 
tribution  (both  polarizations). 


Fig.  29.  Chromaticity  curves  for  the  Airy  and  the  Mie  primary 
fogbows  differ  by  A (u' ,  v')  =  0.00243  (both  polarizations). 


10-|xm-radius  bows  shown  above.  In  Fig.  27,  Crms 
for  the  primary  Mie  and  Airy  fogbows  has  been  re¬ 
duced  43%  from  its  Fig.  7  value.  Figure  28’s  second¬ 
ary  fogbows  show  smaller  improvement;  there  Crms 
has  dropped  —9%  compared  with  that  of  Fig.  21. 
Luminance  ripples  have  vanished  from  Fig.  27’s  Mie 
primary  and  are  all  but  absent  in  Fig.  28’s  second¬ 
ary.38 

Mie-Airy  colorimetric  distances  are  affected  even 
more  than  luminances  by  a  fog  polydispersion.  Fig¬ 
ure  29’s  primary  fogbow  A (uf,  v')  of  0.00243  is  nearly 
70%  smaller  than  that  in  Fig.  9  (10-|xm  primary),  and 
Fig.  30’s  secondary  fogbow  A (ur,  v')  of  0.0071  is  55% 
smaller  than  that  in  Fig.  22  (10-p.m  secondary).  The 
Airy  primary  is  now  within  1  JND  of  its  Mie  equiv¬ 
alent,  and  the  Airy  secondary  probably  would  be  in¬ 
distinguishable  as  well.  Figure  29  also  reveals  how 
the  Mie  color  ripple  structure  disappears  in  the  pri- 


Fig.  30.  Chromaticity  curves  for  the  Airy  and  the  Mie  secondary 
fogbows  differ  by  A (u’,  v')  =  0.00710  (both  polarizations). 
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Table  1.  Summary  of  Chromaticity  Distances  A(u',  v’)  and  mean 
contrast  Crms  between  Mie  and  Airy  Theory  Rainbows  and  Cloudbows® 


Drop  Radii 

(M-m) 

Primary 

(1  and  ||) 

Primary 

A («',  v') 
(L  and  ||) 

Secondary 

cm 

(±  and  ||) 

Secondary 
Mu',  v') 
(±  and  II) 

10 

0.1082 

0.0080 

0.2946 

0.01585 

50 

0.1240 

0.00571 

0.2839 

0.01944 

150 

0.1122 

0.01037 

0.2940 

0.03324 

500 

0.1327 

0.01280 

0.3108 

0.04019 

10-1000 

Cv  =  -0.0227 

0.005885 

Cv  =  -0.1421 

0.01034 

monodis¬ 

perse 

1-40  fogbow 

0.06140 

0.00243 

0.2688 

0.00710 

“The  first  four  rows  are  based  on  Figs.  4-13  and  18-25,  and  for 
these  comparisons,  minima  are  in  italic  type  and  maxima  are  in 
boldface  type.  Entries  in  the  row  labeled  10-1000-jxm  monodis- 
perse  are  averages  over  Figs.  14-17.  The  last  row’s  entries  refer 
to  Figs.  27-30. 


mary  fogbow:  the  Mie  chromaticity  curve  straight¬ 
ens  out  along  the  Airy  arc.  In  Fig.  30,  the 
secondary’s  greater  width  makes  it  easier  to  see  more 
angular  detail,  thus  making  the  Mie  curve’s  chroma¬ 
ticity  wiggles  more  obvious  there.  However,  com¬ 
pared  with  the  two  theories’  stark  colorimetric 
differences  in  Fig.  22,  they  agree  much  more  closely 
in  Fig.  30.  Thus  even  at  the  most  demanding  (i.e., 
smallest)  droplet  sizes,  smoothing  over  a  natural 
polydispersion  largely  eliminates  perceptible  differ¬ 
ences  between  the  Mie  and  the  Airy  theory  bows. 
Furthermore,  my  own  digital  image  analyses  of  nat¬ 
ural  cloudbows  reveal  smooth  luminance  and  chro¬ 
maticity  curves  more  akin  to  Airy  than  Mie 
theory.35’39 

7.  Conclusions 

Table  1  summarizes  the  color  and  luminance  differ¬ 
ences  between  Mie  and  Airy  theory  discussed  above. 
One  of  its  most  remarkable  features  is  that,  for  single 
drops,  Airy’s  largest  color  and  contrast  errors  occur  at 
the  biggest  radii  (500  [xm  here).  However,  if  we  in¬ 
terpret  errors  as  angular  displacements  of  primaries, 
secondaries,  and  supernumeraries  from  their  Mie 
theory  positions,  then  Airy  theory  errors  do  indeed 
increase  with  decreasing  drop  size.  For  both  prima¬ 
ries  and  secondaries,  many  of  the  largest  Mie-Airy 
A (u',  v')  occur  in  Alexander’s  dark  band,  where  low 
luminance  usually  renders  them  invisible.  Al¬ 
though  color  differences  are  often  subthreshold,  con¬ 
trasts  are  not,  and  these  are  largest  for  the  Airy 
secondaries.  In  addition,  angular  shifts  from  the  po¬ 
sitions  of  Mie  maxima  and  minima  are  larger  in  the 
Airy  secondaries  than  in  the  primaries. 

However,  from  a  purely  visual  standpoint,  the  Airy 
monodisperse  bows  (Figs.  14  and  16)  look  more  nat¬ 
uralistic40  than  their  Mie  kin  (Figs.  15  and  17),  and 
they  do  so  because  they  lack  some  of  the  latter’s  spec¬ 
tral  detail.  Furthermore,  if  we  include  a  polydisper¬ 
sion  in  our  calculations  (Figs.  27-30),  Airy  and  Mie 
theories  become  perceptually  indistinguishable  even 
at  cloud-drop  sizes.  Thus,  far  from  being  an  out¬ 


dated  irrelevancy,  Airy  theory  shows  itself  to  be  a 
simple,  quantitatively  reliable  model  of  the  natural 
rainbow’s  colors  and  luminances.  Provided  that  we 
restrict  ourselves  to  spectrally  integrated  luminances 
(or  radiances),  it  can  be  used  to  predict  accurately  the 
visual  appearance  of  most  naturally  occurring  water- 
drop  bows  and  supernumeraries.  Equally  impor¬ 
tant,  we  see  that  Mie  theory  for  monodispersions 
should  not  always  be  the  model  of  first  resort,  as  it 
produces  details  not  seen  in  the  natural  rainbow  and 
does  so  with  much  more  computational  effort  than 
Airy  theory. 
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Colorimetrically  calibrated  versions  of  Figs.  14- 
17  can  be  viewed  at  http://www.nadn.navy.mil/ 
Oceanography/papers/RLee_papers.html.  This  ad¬ 
dress  is  case  sensitive.  If  this  link  is  outdated,  write 
to  the  author  by  post  or  email  (raylee@nadn.navy. 
mil)  for  the  figures’  current  location. 
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